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Murphy-Topel Covariance Matrix for Eviews 


George S. Ford, Chief Economist, The Phoenix Center for Advanced Legal & Economic Public 
Policy Studies, george.ford@phoenix-center.org. 


Abstract. An Eviews program is provided that computes the Murphy-Topel 
Covariance Matrix using Eviews. A dataset is provided to confirm the output, 
which is compared to the output of the procedure using STATA. 


I. Introduction 


The program computes the Murphy-Topel covariance matrix for two-stage estimation 
procedures. Programming is based on the work of Hole (2006). The output created by the 
program matches the output of STATA 8.0 using the data from Greene (2000). I cannot 
guarantee this program will work with all datasets. Both the logit and probit first-stage 
estimates are indicated. 


II. Program 


‘CREATE NEEDED SERIES 

series zhat=na 

series const=1 

‘CREATE X GROUPS 

group x1 age income ownrent selfemp const 

group x2 age income expense zhat const 

‘ESTIMATE FIRST STAGE EQUATION 

equation eq1.binary(d=l) z x1 

‘equation eq1.binary(d=p) z x1 ‘Use this if first stage is probit 

eq1.fit zhat 

“eq1.fit(d) xb “Use this if first stage is probit 

sym v1 = eq1.@cov 

eq1.makegrads 

Ik1 = eq1.Oncoef 

series s1 = grad0!k1 "This matches the 'score' command in Hole (2007) 
delete grad* 

'ESTIMATE SECOND STAGE EQUATION 

equation eq2.count(d=p) y x2 

sym v2 = eq2.@cov 

eq2.makegrads 

Ik2 = eq2.@ncoef 

series s2 = grad0!k2 

scalar zz = eq2.@coef(4) ‘Make sure this gets the coef off the zhat 
delete grad* 

‘CALCULATE C USING SCORES 

series iwpre=s2*s2*zhat*(1-zhat)*zz ‘Eviews does not have 'matrix accum' command 
‘series iwpre=s2*s2*Odnorm(xb)*zz “Use this if first stage is probit 
sym iw = @makediagonal(iwpre) 

group x3 x1 x2 

stom(x1,xx1) 

stom(x2,xx2) 

stom(x3,xx3) 

‘CALCULATE C 

sym bigc = @transpose(xx3)*iw*xx3 

‘CALCULATE R 

series iwpre=s2*s1 

sym iw = @makediagonal(iwpre) 

sym bigR = @transpose(xx3)*iw*xx3 
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'GET DESIRED PARTITION 

matrix smc=@subextract(bigc,6,1,10,5) 

matrix smr=@subextract(bigr,6,1,10,5) 

'FINAL COVARIANCE MATRIX 

matrix m = v2 + (v2*(smc*v1*@transpose(smc) - smr*v1*@transpose(smc) - smc*v1*@transpose(smr)) * v2) 
vector v_se=@getmaindiagonal(m) 

v_se=@sart(v_se) 


III. Data and Output 


Running the program produced the following results. The logit-first stage is presented first, 
then the probit first-stage. The outcomes match those found in Hole (2007). 


First-Stage Logit; Uncorrected Standard Errors. 


Dependent Variable: Y 

Method: ML/QML - Poisson Count (Quadratic hill climbing) 
Date: 11/14/08 Time: 10:08 

Sample: 1 100 

Included observations: 100 

Convergence achieved after 6 iterations 

Covariance matrix computed using second derivatives 


Variable Coefficient Std. Error z-Statistic Prob. 
AGE 0.073106 0.054246 1.347681 0.1778 
INCOME 0.045234 0.174111 0.259797 0.7950 
EXPENSE -0.006897 0.002020 -3.414291 0.0006 
ZHAT 4.632355 3.661768 1.265060 0.2058 
CONST -6.319948 3.930761 -1.607818 0.1079 
Mean dependent var 0.360000 S.D. dependent var 1.010250 
S.E. of regression 0.971226 Akaike info criterion 1.666620 
Sum squared resid 89.61167 Schwarz criterion 1.796878 
Log likelihood -78.33099 Hannan-Quinn criter. 1.719338 


Avg. log likelinood -0.783310 


Corrected Standard Errors (in v_se) 


AGE 0.109629 
INCOME 0.437539 
EXPENSE 0.004265 
ZHAT 10.82668 
CONST 9.661544 
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First-Stage Probit; Uncorrected Errors 


Dependent Variable: Y 

Method: ML/QML - Poisson Count (Quadratic hill climbing) 
Date: 11/14/08 Time: 10:15 

Sample: 1 100 

Included observations: 100 

Convergence achieved after 6 iterations 

Covariance matrix computed using second derivatives 


Variable Coefficient Std. Error z-Statistic Prob. 
AGE 0.080301 0.059369 1.352568 0.1762 
INCOME 0.039716 0.176684 0.224784 0.8221 
EXPENSE -0.006886 0.002013 -3.420316 0.0006 
ZHAT 5.393431 4.229164 1.275295 0.2022 
CONST -7.094362 4.509666 -1.573146 0.1157 
Mean dependent var 0.360000 S.D. dependent var 1.010250 
S.E. of regression 0.970285 Akaike info criterion 1.663875 
Sum squared resid 89.43800 Schwarz criterion 1.794134 
Log likelihood -78.19377 Hannan-Quinn criter. 1.716593 


Avg. log likelinood -0.781938 


Corrected Standard Errors (in v_se) 


AGE 0.150957 
INCOME 0.522169 
EXPENSE 0.00471 
ZHAT 14.91045 
CONST 13.68201 


The dataset used is from Greene (2000), available at 
http://www.stern.nyu.edu/~wereene/Text/Edition6/TableF8-1.txt. 


y Z age income expense ownrent selfemp 
0 1 38 4.52 124.98 1 0 
0 1 33 2.42 9.85 0 0 
0 1 34 4.5 15 1 0 
0 1 31 2.54 137.87 0 0 
0 1 32 9.79 546.5 1 0 
0 1 23 2.5 92 0 0 
0 1 28 3.96 40.83 0 0 
0 1 29 2.37 150.79 1 0 
0 1 37 3.8 777.82 1 0 


0 1 28 3.2 52.58 0 0 
0 1 31 3.95 256.66 1 0 
0 0 42 1.98 0 1 0 
0 0 30 1.73 0 1 0 
0 1 29 2.45 78.87 1 0 
0 1 35 1.91 42.62 1 0 
0 1 41 3.2 335.43 1 0 
0 1 40 4 248.72 1 0 
1 0 30 3 0 1 0 
0 1 40 10 548.03 1 1 
3 0 46 3.4 0 0 0 
0 1 35 2.35 43.34 1 0 
1 0 25 1.88 0 0 0 
0 1 34 2 218.52 1 0 
1 1 36 4 170.64 0 0 
0 1 43 5.14 37.58 1 0 
0 1 30 4.51 502.2 0 0 
0 0 22 3.84 0 0 1 
0 1 22 1.5 73.18 0 0 
0 0 34 2.5 0 1 0 
0 1 40 5.5 1532.77 1 0 
0 1 22 2.03 42.69 0 0 
1 1 29 3.2 417.83 0 0 
1 0 25 3.15 0 1 0 
0 1 21 2.47 552.72 1 0 
0 1 24 3 222.54 0 0 
0 1 43 3.54 541.3 1 0 
0 0 43 2.28 0 0 0 
0 1 37 5.7 568.77 1 0 
0 1 27 3.5 344.47 0 0 
0 1 28 4.6 405.35 1 0 
0 1 26 3 310.94 1 0 
0 1 23 2.59 53.65 0 0 
0 1 30 1.51 63.92 0 0 
0 1 30 1.85 165.85 0 0 
0 1 38 2.6 9.58 0 0 
0 0 28 1.8 0 0 1 
0 1 36 2 319.49 0 0 
0 0 38 3.26 0 0 0 
0 1 26 2.35 83.08 0 0 
0 1 28 7 644.83 1 0 


0 0 50 3.6 0 0 0 
0 1 24 2 93.2 0 0 
0 1 21 1.7 105.04 0 0 
0 1 24 2.8 34.13 0 0 
0 1 26 2.4 41.19 0 0 
1 1 33 3 169.89 0 0 
0 1 34 4.8 1898.03 0 0 
0 1 33 3.18 810.39 0 0 
0 0 45 1.8 0 0 0 
0 1 21 1.5 32.78 0 0 
2 1 25 3 95.8 0 0 
0 1 27 2.28 27.78 0 0 
0 1 26 2.8 215.07 0 0 
0 1 22 2.7 79.51 0 0 
3 0 27 4.9 0 1 0 
0 0 26 2.5 0 0 1 
0 1 41 6 306.03 0 1 
0 1 42 3.9 104.54 0 0 
0 0 22 5.1 0 0 0 
0 1 25 3.07 642.47 0 0 
0 1 31 2.46 308.05 1 0 
0 1 27 2 186.35 0 0 
0 1 33 3.25 56.15 0 0 
0 1 37 2.72 129.37 0 0 
0 1 27 2.2 93.11 0 0 
1 0 24 4.1 0 0 0 
0 1 24 3.75 292.66 0 0 
0 1 25 2.88 98.46 0 0 
0 1 36 3.05 258.55 0 0 
0 1 33 2.55 101.68 0 0 
0 0 33 4 0 0 0 
1 1 55 2.64 65.25 1 0 
0 1 20 1.65 108.61 0 0 
0 1 29 2.4 49.56 0 0 
3 0 40 3.71 0 0 0 
0 1 41 7.24 235.57 1 0 
0 0 41 4.39 0 1 0 
0 0 35 3.3 0 1 0 
0 0 24 2.3 0 0 0 
1 0 54 4.18 0 0 0 
2 0 34 2.49 0 0 0 


0 0 45 2.81 0 1 0 
0 1 43 2.4 68.38 0 0 
4 0 35 1.5 0 0 0 
2 0 36 8.4 0 0 0 
0 1 22 1.56 0 0 0 
1 1 33 6 474.15 1 0 
1 1 25 3.6 234.05 0 0 
0 1 26 5 451.2 1 0 
0 1 46 5.5 251.52 1 0 
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